Not Much Helicity is Needed to Drive Large Scale Dynamos 
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Understanding the in situ amplification of large scale magnetic fields in turbulent astrophysical 
rotators has been a core subject of dynamo theory. When turbulent velocities are helical, large scale 
dynamos that substantially amplify fields on scales that exceed the turbulent forcing scale arise, but 
the minimum sufficient fractional kinetic helicity fh,c has not been previously well quantified. Using 
direct numerical simulations for a simple helical dynamo, we show that fh,c decreases as the ratio of 
forcing to large scale wave numbers ftF /k m i n increases. From the condition that a large scale helical 
dynamo must overcome the backreaction from any non-helical field on the large scales, we develop 
a theory that can explain the simulations. For Uf /kmin > 8 we find fh,c IS 3%, implying that very 
small helicity fractions strongly influence magnetic spectra for even moderate scale separation. 



>> 
q 

O 
Of 



>: 

OV 

m 

O 

rn 

O 



Introduction The origin of magnetic fields in turbulent 
astrophysical rotators such as stars, galaxies [l[ and ac- 
cretion disks has been a long standing topic of research. 
A particular challenge has been to understand the origin 
of fields on scales that are large compared to those of any 
underlying turbulence [J|3|- 

That the large scale field of the sun reverses every 11 
years reveals that such stellar fields cannot be simply 
the residual of flux freezing from the primordial mate- 
rial and must be amplified in situ. Complcmcntarily, the 
continuous processing by supcrnovac driven turbulence 
in galaxies likely renders the role of any primordial fields 
to be simply seed fields whose in situ processing must be 
understood to account for the observed present day large 
scale fields in galaxies. The presence of astrophysical 
jets from accretion engines also highlights the presence 
of large scale fields in accretion disks, and accretion disk 
simulations Q commonly show the in situ generation of 
large scale magnetic fields that reverse on cycle periods 
of tens of orbit times. 

The study of in situ field amplification in the pres- 
ence of velocity flows is the enterprise of dynamo theory. 
Small scale dynamos (SSDs) , in which turbulent velocity 
flows amplify fields at or below scales of the forcing [1,0], 
can be distinguished from large scale dynamos (LSDs) in 
which magnetic fields arc amplified on spatial or temporal 
variation scales larger than the scales of the underlying 
forcing. LSDs and SSDs are often contemporaneous and 
interactive (see e.g. (3,01) but LSDs arise only when tur- 
bulent velocities are sufficiently helical [H, Q . There has 
been little previous work, however, on determining the 
minimum sufficient helicity to incite LSD action and this 
is the topic of the present paper. Astrophysical flows 
are unlikely to be 100% helical in environments where 
LSDs are presumed; the galaxy for example is estimated 
to have helicity of < 10%. Thus the basic question of 
how much helicity is required in even the simplest LSDs 
is important in assessing the potential ubiquity of LSDs. 



The standard 20th century textbook [T(| kinematic 
approach to LSD theory has been classical mean field 
(MFT) which features the a-effect: 7 = \a\k - f3k 2 , 
where 7 is the exponential growth rate in the kinematic 
dynamo regime (presuming that any Lorentz-forcc feed- 
back is negligible), k is the wavenumber of magnetic field 
growth, the a— effect is proportional to kinetic helicity 
H v — (v • u>), with v the velocity and uj — V x v the 



vorticity, and /3 is the turbulent eddy diffusivity 1C | 
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Such mean-field theory has been used to model solar 
stellar, and galactic observations, as well as laboratory 
plasma dynamos [l2| , and Geo-dynamos [l3[ . 

But the kinematic approach to LSD theory is incom- 
plete. Although many astrophysical rotators have differ- 
ential rotation and open boundaries, substantial progress 
in going beyond the kinematic theory has emerged from 
studies of the closed volume "a 2 " helical dynamo with- 
out shear, in which the evolution of an initially weak 
seed field is subject to helical velocity forcing. The 
a 2 dynamo was first tackled semi-analytically H us- 
ing a spectral integro-differential model with an Eddy 
Damped Quasi-Normal Markovian (EDQNM) closure, 
consistently tracking the magnetic helicity. It was shown 
that the actual driver of large scale magnetic field growth 
is not just the kinetic helicity, but the residual helicity, 
Hb = H v — Hj where the current helicity Hj = (j • b) 
and j = V x b is the current density. 



This a 2 dynamo in a periodic box was simulated [14j by 
forcing with kinetic helicity at wavenumber hp = 5k m i n 
(k m i n = 1 was the smallest wavenumber of the flow) . The 
large-scale (k < 5) field grew as expected from Ref. M. 
Subsequently, a two-scale a 2 LSD was developed [Fa ]; 
it incorporated magnetic helicity evolution using a sim- 
pler closure than EDQNM and showed even a two-scale 
nonlinear theory predicts the evolution and saturation of 
LSD growth observed in [l4| . Driving with kinetic helic- 
ity initially produces a large scale helical magnetic field, 
but the near conservation of magnetic helicity leads to a 
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compensating small scale magnetic (and current) hclic- 
ity of opposite sign. This counteracts the kinetic helicity 
driving in the large scale field growth coefficient, and 
quenches the LSD, as proposed in Q. 

There has also been a plethora of work on the SSD. In a 
periodic box with a weak initial seed field and non-helical 
forcing, the stochastic line stretching produces negligible 
field growth above the forcing scale. Simulations of non- 
helical SSDs without large scale shear show that the total 
magnetic energy is amplified to near equipartition with 
the total kinetic energy not only in the kinematic regime 
as predicted by but also in the saturated regime for 
large magnetic Prandtl number Pm = v/t), where v and 
rj are the viscosity and magnetic diffusivity @, HH, E3 , 
as reviewed in pj. Astrophysical plasmas such as the 
Galactic interstellar medium do not seem to exhibit this 
pile-up [lj]. 

To address this disparity between simulations of non- 
helical SSDs and how conditions favorable for LSD 
growth might influence the magnetic spectrum on both 
large and small scales, results for dynamos in a periodic 
box forced with different amounts of fractional kinetic he- 
licity fh (a dimensionlcss measure of the degree of align- 
ment between the velocity and the vorticity of the forcing 
function), were studied [U^. It was found that the mag- 
netic spectrum above and below the forcing scale were 
contemporaneously affected by a sufficient fh- The large 
scale field grew, and the magnetic spectrum at large wave 
numbers steepened. For fh = l and fh = 0, the results 
of Q an d [13] were respectively recovered. 

But the restriction in to a forcing scale of kp = 5 
and resolution of 64 3 grid points left key unexplored ques- 
tions. In particular, the minimum fh for LSD action, 
fh,c, could not be determined as a function of kp . The 
smaller this minimum, the potentially more ubiquitous 
LSD conditions are in astrophysics. Here we perform 
much higher resolution simulations for fractionally heli- 
cal dynamos and quantify how fh.c depends on kF/k m i n . 
We also develop a theory that correctly predicts the de- 
pendence, seen in the simulations. 

Equations and set-up- The incompressible MHD equa- 
tions for velocity v and magnetic field b are: 

d,y + w x v = j x b - Vp + zA7 2 v + F 
d t A = v x b - + 77V 2 A 

V-v = 0, V-A = 0. (1) 

The total pressure divided by the constant (unit) density 
p and the potential <j) are obtained self-consistently to 
ensure incomprcssibility and the Coulomb gauge. The 
Reynolds number is Re = U rms Lo/v, with U rms and 
L a = 2tt J E(k)k~ 1 dk/ J E(k)dk the r.m.s. velocity and 
the integral scale respectively; the magnetic Reynolds 
number is defined as Rm = U rms Lo/ri. In the follow- 
ing, E denotes the total energy, and E v and Ef, denote 
the kinetic and magnetic energy respectively. 

We employ a well-tested pseudo-spectral code that uses 
a hybrid parallelization, combining Message Passing In- 
terface (MPI) and OpenMP [20|]. The computational box 



has size [27r] 3 , and wave numbers vary from k m i n = I to 
kmax = N/3 using a standard 2/3 de-aliasing rule, where 
N is the number of grid points per direction. 
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FIG. 1. (Color online) Re-dimensionalized magnetic (solid) 
and kinetic (dashed) energy spectra after 90r for run 4-60 
(thick) black forced at kp =4 and for run 3-60 (thin red/light 
gray) with kp = 3. At a fixed fh — 60% here, increased scale 
separation provides for the transition between SSD and LSD. 

The forcing applied at kp is F = F# + cFa', Fa is 
an ABC flow at kp, and F# is the sum of all harmonic 
modes with k = kp and random phases. We choose c 
for a given fractional helicity of F, \fh\ < 1, with fh = 
(F-a; F )((|F| 2 )(|u; F | 2 ))- 1 /2 ! where w F = V x F. The en- 
tire forcing has random phases applied, with a correlation 
time t cor = 0.1. In practice, the ratio of helical to non- 
helical forcing magnitudes is c ~ [fh/0- ~ fh)] 1 ^ 2 = Rh- 
The kinetic helicity at fc F is typically within 25% of 
fh- By choosing dimensional length and time constants 
lo oc kp and to (fixed), varying kp in our simulations cor- 
responds to dimensionalizcd physical systems described 
by Eqs. ([T]), where the forcing scale is constant and the 
system size increases oc kp. The dimensionless velocity 
v rms and forcing are oc kp 1 , and the diffusivity oc kp 2 . 

A hydrodynamic state is evolved for five forcing-scale 
eddy turnover times, r, before a magnetic seed field at 
k = k seec i is introduced. In the hydro steady state, the 
resulting t = 2~K\kpU rms \ « 4.2. In all simulations, 
dimensionalized viscosity is constant, k F v = 2.412-10~ 2 , 
and, arbitrarily, P^j = 4 so that P^j > 1 while limiting 
the computational cost (see Table |l] for further details). 

Simulation Results- Table [J summarizes our runs. The 
kinetic and magnetic energy spectra after t = 90r are 
displayed in Fig. [I] for runs 4-60 (with kp = 4) and 3-60 
(with kp = 3). Both have //j=60%, and for both, the 
small-scale fields grow; only for kp = 4 does k = 1 grow. 

Figure [2] shows the growth of magnetic energy in the 
k = 1, k = 6, and total over all modes for run 3-80. The 
evolution exhibits an early phase in which both modes 
grow at the same rapid rate, with jssd ~ 0.33 ~ t _1 , 
followed by a slow growth of the k = 1 mode and a sat- 
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TABLE I. Parameters: Runs are labeled by the forcing wavenumber Icf followed by the percentage of helicity in the forcing fh\ 
Rm and fc seet j are defined in the text (Re = Rm/4)- The SSD growth rate is "/ssd- Ef, is the magnetic energy and Hb is the 
magnetic helicity, both at later times, between 60r and 130r; the "/" is for fluctuating, and "g" is for growing exponentially 
(indicating a helical dynamo). Forcing wave numbers are fcf = 2, 3, 4, 5, 6 and 8 for runs on grids of 192 3 , 256 3 , 384 3 , 432 3 , 512 3 
and 768 3 points respectively. *Run 3-80 was pursued until t ~ 300r, at which time Eb ~ 0.3 and Hb ~ 10. See Fig. [2] 
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FIG. 2. Magnetic energy density Eb(k) for k = 1 (solid line), 
k = 6 (dashed) and total (dotted) versus time for run 3-80. 
The gray line indicates the fit 7 = (8.6 ± 1.4) 10~ 3 to the solid 
curve for growth of the k — 1 mode after SSD saturation. 



FIG. 3. Growth rate 7^=1 of Eb(k — 1) versus fractional 
helicity. From Table 0] x (pink, kp = 2), + ( green, kF = 3), 
□ (cyan, kF = 4), A (blue, = 5), o (red, kF = 6), grey, * 
(kF = 8) and least-squares linear fits (dashed). 



uration for the k = 6 mode. The k = 1 mode accounts 
for nearly 10% of Eb by 100Y and the growth rate slows, 
but has not fully saturated by 300r. The growth rate of 
magnetic energy at k = 1 during the SSD phase is nearly 
the same for all of our runs, and is insensitive to fh, Rm, 
and kp. Sensitivity to fh emerges once the SSD regime 
ends. 

The 7fc=i growth rates (for Eb(k = 1)) that immedi- 
ately follow the SSD phase (see TablcHJ are shown in Fig. 
[21 This LSD growth regime occurs only when fh > fh,C] 
the LSD growth rate varies linearly with fh for a fixed 
kp/kmin {pt oc H v oc fh). Least-squares fits are dashed 
lines in Fig. [3] (the y-intercept, (3 ~ (k min /kp) 2 - 2 ). The 
short exponential growth phase of 4-80 makes for an in- 
accurate measure of 7fc=i; it is thus excluded from the 



fit. As kF/k m i n increases, fh,c decreases. 

The LSD exponential growth of Eb(k = 1) for fh > 
fh,c is accompanied by a k = 1 growth of magnetic helic- 
ity. Studies of the k = 1 growth for fh = 1 in a two-scale 
approach [3, [l5| for a Hr driven dynamo [j| suggest two 
phases of k m i n = 1 mode growth after the SSD regime: 
one phase that is largely independent of Rm , and a sub- 
sequent Rm dependent asymptotic regime. The former 
phase has growth consistent with our 7^=1 phase. In all 
these runs, fh,c decreases with increasing kp/k m i n , as 
displayed in Fig. 01 For the largest kp/k m i n (=8) case, 
fh ~ 3% is sufficient. In Fig. |4] the error bars and 
fh,c arc calculated as follows: The x-intercept from the 
least-squares linear fits, "fk=i = mfh + b shown in Fig. [31 
determine our estimate of fh,c = —b/m for Fig. 01 The 
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FIG. 4. fh,c from least-squares fits versus kF/k m i n (symbols 
as in Fig. [3J). Dashed line is best fit to Eq. giving C — 
0.21 and £ = 0.46. Dotted line is kinematic MFT prediction 
fh.c = /3k/(\a \k F ). 



I— a uncertainties for b and m are then propagated for 
fh.G ■ Vf h , c = b/m^(o- b /b) 2 + (cr m /m) 2 . 

Theoretical prediction for fh,c~ The following predic- 
tion for fh,c is based on the principle that LSD helical 
field growth at k = k m i n beyond the SSD phase requires 
helical velocity forcing to overcome the Lorentz force at 
k = k m i n at the end of the kinematic SSD phase. For 
that time, we assume the magnetic energy at k m i n < kp 
to be B^ in ~ B F {k min /kp)^ , where Bp is the magnetic 
field at kp , and ^ — 1 is the slope of the magnetic energy 
spectrum on a log-log plot. The associated Lorentz force 
is then M n h(fh)B F kp(k mirl /kp)^ +1 , where the function 
M n h(fh) < 1 accounts for the contribution from only 
non-helical magnetic energy. 

The available helical velocity forcing that must over- 
come this Lorentz force is only a fraction of the helical 
forcing at k = kp: At early times when magnetic helicity 
is nearly conserved, the forcing not only sources magnetic 
helicity at k = k m i n but also an oppositely signed, equal 
in magnitude, magnetic helicity at k = k ss > kp. The as- 
sociated ratio of helical magnetic energy growth at fc m j n 
to that at k ss is then ~ k min /k ss < 1. The helical force 
that needs to exceed the Lorentz force at k m i n to ini- 
tiate growth is thus ~ K h (f h )v 2 ms kp(k min /k ss ), where 
the function Kh(fh) < 1 accounts for only kinetic helical 
forcing. 

Balancing the aforementioned forces assuming 
Kh/M n h = Rh(fh) (which is consistent with our data), 
and assuming fh = fh.c, then gives 



fh,C 



1 



l + C 2 (k F /k min ) 



2?+2 ■ 



where C = (k mm v 2 ms )/{k ss B F ) 



i/k s 



(2) 



Fig- 



ure |3] shows the data and the best fit using Eq. ©; 
£ = 0.46 ~ 1/2 is found. This yields the prediction, 
fh,C ~ (k F /k mm )~ 3 as k F /k min -> oo. Note that taking 



the limit of infinite scale separation, we have a LSD with 
zero helicity (but only fluctuations, as in [2l|). 

The theory above, which considers the Lorentz force 
backreaction from the large scale field, can be contrasted 
with the prediction from the purely kinematic theory 
of the standard a 2 dynamo which does not include any 
Lorentz forces. Using the formula presented in the intro- 
duction for the growth rate at k m i n , and the definition of 
fh, the critical fractional helicity for the kinematic the- 
ory would be f h ,c = /3k min / (\a \k f ) where a = a/ f h . 
This formula is shown as the dotted line in Fig. 2] and 
does not fit the data very well, highlighting the impor- 
tance of including the Lorentz force. This does not imply 
the kinematic theory is irrelevant however. For values of 
fh >> fh.c the kinematic theory should be applicable to 
estimating the early time growth rate because the driv- 
ing helicity overwhelms the backreaction associated with 
the weak large scale field produced by the SSD in that 
regime. 

Discussion of LSD growth and saturation- At large Rm, 
SSD action produces field at all scales, potentially pre- 
cluding a scale separation between the mean magnetic 
field and velocity fluctuations (an essential assumption 
to derive a 2 MFT) [22| ■ The SSD magnetic energy spec- 
trum at scales above the forcing scale produces less mag- 
netic energy the larger the scale [f|. At large enough 
scales, the magnetic energy production from the SSD will 
be negligible, and scale separation becomes a meaning- 
ful concept (see thin red/light gray, solid line in Fig. [TJ. 
This helps justify the mean field approach to LSDs. 

In the mean field, two-scale approach, once the small 
scale magnetic helicity has grown as a result of magnetic 
helicity conservation to be large enough such that the 
associated small scale current helicity backreacts on the 
driving kinetic helicity, the a 2 dynamo eventually slows 
to Rm —dependent growth rates and ultimately saturates 
completely Previous studies have typically focused on 
the fh = 1 case [HI, HH. In this paper, we have not run 
enough simulations long enough to determine how strong 
the large scale field gets before its evolution reaches the 
Rm dependent regime. However, if cases with fractional 
helicity fh.c < fh < 1 saturate by direct analogy to the 
fh = 1 cases studied in previous work, then the value 
of the large scale magnetic energy reached just before 
the Rm dependent regime emerges would be expected to 
be simply proportional to an extra factor of fh, namely 

2 



B 2 



fh{ki/kf)(U 2 ms ). Similarly, for asymptotically sat- 



urated steady state at very late times, we would expect 

B ~ fh(kf / k\)(U 2 ms ) . Note that in the fh = 1 case, 
the latter similarity highlights the fact that that supcr- 
cquipartition field strengths (with respect to the total 
kinetic energy) are able to grow by the end of the non- 
linear, saturated regime for fully helical a 2 dynamo. 

Note however that the Rm dependent regimes of the 
a 2 dynamo are largely irrelevant for astrophysical objects 
which have such large Rm that something else probably 
happens before these regimes are reached. Open bound- 
aries and helicity fluxes are ingredients that have to be 
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considered in realistic systems. In addition, real astro- 
physical dynamos have large scale shear, which amplifies 
the total large scale field beyond its purely helical value. 
More work is needed to determine the strength of the 
large scale fields produced by fractionally helical LSDs. 

Conclusion- Only a minuscule amount of fractional he- 
licity is required for LSD action at even modest astro- 
physically relevant scale separations. For fh > fh,Ci the 
k = 1 field grows and the small scale spectrum steepens 
(see Fig. [T] and [IH). This may be important because 
our result that fh.c i5 3% for k F /k min > 8, offers a basic 
principle for potentially reconciling a disparity between 



a pile up of small scale magnetic energy in large Pm non- 
helical dynamo simulations [f| and Galactic observations 
[lH . Our results also suggest that large scale separations 
should be a priority in designing laboratory experiments 
to measure LSD action [2j. 
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